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An aerothermodynamic analysis of the forebody aeroshell of the Stardust Sample Return Capsule is carried out 
by using the axisymmetric viscous shock-layer equations with and without fully coupled radiation and ablation. 
Formulation of the viscous shock-layer equations with shoulder radius as the length scale and implementation 
of the Vigneron pressure condition allow resolution of the flowfield over the shoulder. With a predominantly 
supersonic outflow over the shoulder, a globally iterated solution of viscous shock-layer equations can be obtained. 
The stagnation-point results are obtained along a specified trajectory, whereas detailed calculations along the body 
are provided at the peak-heating point. The equilibrium calculations with ablation injection are the focus of the 
present study because of the lack of a general chemical nonequilibrium analysis that accounts for both surface 
and flowfield effect. The equilibrium calculations also provide a simple way to conserve surface (and flowfield) 
elemental composition for the current small ablation injection rates, where the surface elemental composition is 
a mixture of freestream and ablator elements. Therefore, the coupled laminar and turbulent flow solutions with 
radiation and ablation are obtained by using the equilibrium flow chemistry, whereas a nonequilibrium chemistry 
model is used for solutions without ablation and turbulence. Various computed results are compared with those 
obtained by the other researchers. 
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Nomenclature 

= mass fraction of species i 

- mass fraction of element k 

- frozen specific heat of mixture, J/kg-K 
= binary diffusion coefficient, m 2 /s 

= enthalpy of undecomposed ablation 
material, J/kg 

= enthalpy of species /, J/kg 
= enthalpy of species /, h\j 
= thermal conductivity of mixture, W/m-K 
= thermal conductivity of mixture, K/p n iC P ,oo 
= Lewis number, pD^Cp/K 

- mass injection rate, kg/m 2 -s 

= coordinate measured normal to body, m 

- coordinate measured normal to body, n/Rc i 
= shock standoff distance, m 

= Prandtl number, pC P /K 
= pressure, N/m 2 
= heat flux, W/m 2 
= heat flux, q/PooV^ 

- q cond <?conv “F *7diff> W/m 

= net radiative heat flux in n direction, q { + } - q { ~\ 
W/m 2 

= component of radiative flux toward shock, W/m 2 
= component of radiative flux toward wall, W/m 2 
= q c + q r > W/m 2 
= shoulder radius, m 
= nose radius, m 

- coordinate measured along surface, m 
= temperature, K 

= temperature, T / T n f 
= radiative equilibrium wall temperature, K 
= reference temperature, V£jC Pf<xn K 
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sublimation temperature, K 
time, s 

freestream velocity, m/s 
curve-fit coefficients for sublimation 
temperature, T suh 

curve-fit coefficients for heat of ablation, A H a 
heat of ablation, MJ/kg 
char emmissivity, 0.9 
Reynolds-number parameter, 

(Pmt/poo Voo#ci ) 1/2 
viscosity of mixture, N-s/m 2 
viscosity of mixture, 
reference viscosity, p(T n f), N-s/m 2 
density of mixture, kg/m 3 
Stefan-Boltzmann constant, 

5.668 x 1CT 8 W/(m 2 -K 4 ) 

Subscripts 

A = ablator 

abl = ablation 

cond = conduction 

conv = convection 

diff = diffusion 

Eq = equilibrium 

i - ith species 

j = yth species 

k = Jkth element 

r = radiation 

w = wall value 

- = values for solid ablation material at surface 

oo = freestream value 

Introduction 

T HE Stardust mission, 1 part of NASA’s Discovery Program, 
plans to fly a spacecraft through the tail of the comet Wild-2 
and bring samples of cometary material as well as interstellar dust 
to Earth for analysis. The collected cometaiy particles and the dust 
will be contained in the Stardust Sample Return Capsule (SRC), 
which must survive an intense Earth entry heating. At 12.6 km/s the 
SRC entry is the fastest ever attempted into the Earth’s atmosphere. 
This paper focuses on the aerothermodynamic issues concerning the 
flow environment around the SRC forebody during such an entry. 

Since Stardust was scheduled for launch in early 1999, the work 
presented here was not available in time to impact the SRC design. 
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However, a better understanding of the SRC entry environment and 
the computational tools employed for its analysis will help in pro- 
moting a better and more efficient design of the thermal protection 
system (TPS) for future sample return vehicles such as MUSE S-C, 
Genesis, Champollion (DS-4), and Mars Sample Return. Improved 
computational tools will also be useful in postflight evaluation of 
the TPS and other measured quantities. 

Previous Work 

Recently, calculations have been done 2 for the Stardust SRC fore- 
body TPS with an axisymmetric Navier-Stokes flow solver, loosely 
coupled to the radiation and ablation modules. The thermochem- 
ical nonequilibrium flowfield calculations with ablation are based 
on an 18-species chemical model. The ablation boundary condition 
(namely, the blowing rate, species mass fractions, and wall temper- 
ature) for the flowfield solution are obtained iteratively by assuming 
the surface composition to be in equilibrium at the temperatures and 
pressures predicted from a material response code (with inputs of 
wall heat transfer rate and pressure from the flowfield solution). 

Employing a methodology similar to that of Ref. 2, Ref. 3 has 
recently obtained stagnation-point heat transfer rates for the Pioneer- 
Venus probes, 4 whose flight environment resembles that of current 
sample return vehicles. 2 ’ 5 * 6 A 17-species nonequilibrium chemistry 
model is used for the shock-layer flow. The pyrolysis gas compo- 
sition at the surface is obtained by assuming that the surface is in 
equilibrium at the local temperature and pressure. 

A recent analysis 7 of the MUSES-C asteroid sample return mis- 
sion has considered a 19-species nonequilibrium chemistry model 
(with the thermal equilibrium assumption) both for the shock-layer 
flow and the ablator surface for a true ablation calculation. The 
chemistry model consists of 1 1 air species and 8 carbon-containing 
species. The hydrocarbon species are not included, and the pyroly- 
sis process is not considered to keep the analysis simple. Because 
the 19-species finite rate chemistry model is implemented both at 
the surface and in the shock layer, it includes all of the species con- 
sidered throughout the computational domain. Thus, the analysis is 
consistent both at the surface and through the shock layer for the 
chemistry model considered. 


condition. Only at lower temperatures (i.e., T < 2000 K) would a 
fully catalytic boundaiy condition (with complete recombination to 
freestream value) be realistic to use. 

Analysis 

Flowfield Model 

The viscous shock-layer (VSL) equations employed are those of 
a multicomponent reacting-gas mixture under conditions of chem- 
ical nonequilibrium 9 and equilibrium 8 10 with thermal equilibrium. 
These equations are the same as those given in Refs. 9 and 10, and, 
therefore, they are not given here. The chemistry model, boundary 
conditions, and thermodynamic and transport properties employed 
are similar to those of Refs. 9-1 1, whereas the ablation injection, 
radiative transport, and turbulence models (used with equilibrium 
chemistry only) are those of Refs. 8 and 10-13. These models, 
boundary conditions, and the properties are briefly described here. 

Chemistry Model 

For calculations of airflow over a nonablating surface, an 1 1 -spe- 
cies (N 2 , 0 2 , N, O, NO, NO+, e^ , N + , 0+, N+, and 0+) chemistry 
model is used for nonequilibrium calculations, whereas a 9-species 
(N 2 , 0 2 , N, O, NO, NO + , e _ , N + , and 0 + ) chemical model is used 
for the equilibrium flow. For the equilibrium ablation injection cal- 
culation, 20 chemical species are used: the 7 equilibrium air species 
(without NO and NO + ) plus C, C 2 , C 3 , CO, CN, C 2 H, C 3 H, C 4 H, 
C 2 H 2 , C + , H, H 2 , and HCN. The equilibrium composition is deter- 
mined (for a given temperature, pressure, and elemental composi- 
tion) by using the free-energy minimization method of Ref. 14. 

Boundary Conditions 

The boundary conditions at the shock are obtained by using the 
Rankine-Hugoniot relations. The flow behind the shock is assumed 
to be in chemical equilibrium or frozen at the freestream composi- 
tion ior equilibrium and nonequilibrium calculations, respectively. 9 
No-slip continuum boundary conditions are employed at the sur- 
face. The surface temperature with no ablation injection is assumed 
to be the radiative equilibrium wall value obtained from 


Present Work 

In the present work the Stardust SRC entry flowfield is inves- 
tigated by assuming complete thermal equilibrium. Fully coupled 
radiation solutions with and without ablation injection are obtained 
by using an axisymmetric viscous shock-layer method by assuming 
chemical equilibrium both in the flowfield and at the surface. The 
elemental continuity equations are solved iteratively for each ele- 
ment to determine the appropriate mix of ablative and freestream 
elemental composition at the surface (adjacent to the ablator) and 
through the flowfield. With information concerning the elemental 
composition, pressure, and temperature, the species concentration 
at the surface (adjacent to the ablator) and through the flowfield is 
obtained by using the free energy minimization procedure. 

With the calculation of both the flowfield and surface ablation 
species from the equilibrium assumption, the transition from abla- 
tion to freestream elements (and the corresponding species) is ob- 
tained continuously 8 through the solution of elemental continuity 
equations. Further, the fully equilibrium calculations also bypass the 
entire discussion about governing processes and intermediate steps 
concerning the number of species, reaction mechanisms, and the 
associated reaction rates, especially for the complex flowfields with 
ablation. Therefore, the flowfield calculations with the thermochem- 
ical equilibrium model would appear adequate for the TPS design 
(for the convective-do mi nated heating environment), especially if 
the ablation boundary condition is specified with the equilibrium 
assumption. Although a nonequilibrium overall solution including 
surface and flowfield would be desirable, a full equilibrium calcula- 
tion further provides a conservative estimate of the surface heating 
with the desired safety factor for the TPS design. 

Computed results for an ablating surface include those with the 
laminar flow assumption as well as those for a fully turbulent flow 
immediately downstream of the stagnation line. For a nonablating 
surface results have been obtained with nonequilibrium chemistry 
and fully catalytic 2 as well as equilibrium catalytic wall boundary 



For ablation injection cases steady-state ablation is assumed. How- 
ever, the species surface concentrations, ablation rates, and surface 
temperatures, in general, can be obtained from a material response 
code (such as FIAT of Ref. 2), by employing input surface heat flux 
and pressure from an equilibrium flowfield code. For the surface 
ablation cases considered in the present study, an energy balance at 
the flowfi eld-ablator interface gives the coupled mass injection rate 
for quasi-steady ablation: 


The surface temperature for the present calculations with abla- 
tion injection is that at which the quasi-steady ablation occurs. 
The expression used for surface temperature for the Phenolic 
Impregnated Ceramic Ablator (PICA) 55 has been obtained by curve 
fitting these values (computed from the charring material and ab- 
lation thermal 1617 response code) in the pressure range 0.001 atm 
< p < 1.00 atm. The elemental composition of PICA is similar to 
that of a carbon-phenolic ablator. It is, however, less dense and has 
much lower thermal conductivity. The expressions for the sublima- 
tion temperature and heat of ablation for PICA (with 92% carbon, 
4.9% oxygen, 2.2% hydrogen, and 0.9% nitrogen by mass) are 

5 5 

Taw = ^^ijOogQ)^' 1 + (log p w ) y^a 2J (log C A ) J ~ l 

i 

5 

Oog (log ~ * (3) 

j = i 
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Table 1 Coefficients for ablation temperature for PICA 





j 



Coefficients 

1 

2 

3 

4 

5 

<*U 

«2 J 
«3 

3790.0 

329.94 

20.386 

86.795 

-66.703 

-17.654 

-2980.0 

-1524.6 

-268.62 

-8250.2 

-4340.9 

-771.00 

-7631.7 

-3885.7 

-684.89 


where p w is the wall pressure in atmospheres and C A is the ablator 
mass fraction at the surface. The values for a Uj are given in Table 1 . 
For the case when the gas species adjacent to the surface are caused 
solely by the ablation species (i.e., C A = 1 .0), Eq. (3) with the values 
of a\j in Table 1 gives 

r abl = 3790.0 + 329.94(log p w ) + 20.386(log p w ) 2 (4) 


The surface temperature and the coupled mass injection rate are cal- 
culated by iterating the solution of the governing flowfield equations 
and the boundary conditions. 

For ablation injection the elemental concentrations at the wall 
are governed by convection and diffusion and are obtained from the 
equation 




For the radiative transport calculations the bow shock is considered 
transparent, and the freestream is considered cold and transparent. 
Therefore, the precursor effects are neglected. Further, the body 
surface is assumed to be gray with a reflectivity of 0. 1 , emissivity of 
0.9, and transmissivity of 0. The energy reradiated from the surface 
is included both in the radiation transport calculation as well in 
the surface energy balance [Eq. (2)]. The net radiative flux can be 
represented as 


At the surface 


T 

1 

+ 

II 

(6) 

a — so T ^ 

*ir,w tu 1 w 

(7) 


The heat transferred to the wall because of conduction, diffusion, 
and convection is 


~Pc, w — & Poo 


-3 T _ Le ^ - 3 C f 

3 h Pr ^ dn 


m 

Poo ^00 


NS 

^ \ Cjhj Cf_/i/-)u/ 

i = 1 


where NS is the number of species. 


( 8 ) 


Catalytic Wall Conditions 

For nonablating, nonequilibrium flow the following three cat- 
alytic wall boundary conditions are used. 

1) Noncatalytic wall (NCW): Because no reactions occur at the 
surface in this case, the mass-fraction gradients for all species are 
zero at the surface, i.e., 

m.- 

2) Equilibrium catalytic wall (ECW): The wall catalyzed reac- 
tions are assumed to occur at an infinite rate, and, therefore, the 
species mass fractions at the wall are those corresponding to their 
local equilibrium values, i.e., 

(C i ) w = (C i ) J * l = f(p w ,T w ) (10) 

3) Fully catalytic wall (FCW): The gas species at the surface are 
assumed to recombine to the freestream composition, i.e., 

(C,)* = (C,)« (11) 

At low surface temperatures the surface condition of Eq. (10) re- 
duces to Eq. (11). 



Fig. 1 Geometry of SRC; dimensions in millimeters. 


Radiative Transport 

The radiation transport code RADICAL 1819 has been used to 
compute radiative heat flux q r . This code accounts for the effects of 
nongray self-absorption and includes the molecular band, contin- 
uum, and atomic line transitions. The ultraviolet properties for C3 
are taken from Ref. 20. The radiative transport is fully coupled with 
the flowfield solutions for equilibrium chemistry. The nonequilib- 
rium total heat transfer rate (for a nonablating surface) is obtained by 
adding the equilibrium radiative component to the nonequilibrium 
conductive and diffusive components. 

Thermodynamics and Transport Properties 

Thermodynamic properties for specific heat, enthalpy and free en- 
ergy, and transport properties for viscosity and thermal conductivity 
are required for each species considered. Values of these properties 
are obtained by using polynomial curve fits of Refs. 8 and 1 1. The 
equilibrium composition is determined by a free-energy minimiza- 
tion calculation procedure of Ref. 14. Mixture viscosity is obtained 
by the method of Armaly and Sutton, 21 and mixture thermal con- 
ductivity is computed by the Mason and Saxena 22 relation. Both 
constant and variable Lewis number values are employed in the 
computations. However, the results presented here are for a variable 
Lewis number based on the effective diffusion coefficient D m for 
the mixture. 8 

Turbulence Model 

A two-layer, eddy-viscosity, Cebeci-Smith turbulence mod- 
e j 13, 23,24 j s employed in the present investigation. Reference 24 gives 
a detailed description of the model and various expressions for it. 
The boundary-layer edge definition used in the current study is based 
on an index of diffusion, conduction, and dissipation. 25 The transi- 
tion to turbulent flow is assumed to occur instantaneously at the first 
grid point downstream of the stagnation point. The turbulent Prandtl 
and Lewis numbers are assumed to be 0.9 and 1.0, respectively. 

Solution Procedure 

The method used to solve the nonequilibrium and equilibrium 
VSLequations is a spatial-marching, implicit, finite-difference tech- 
nique, 8 - 13 which includes coupling of the global continuity and nor- 
mal momentum equations and use of the Vigneron pressure condi- 
tion in the subsonic region (which covers a large part of the forebody 
of Stardust capsule shown in Fig. 1). The shoulder radius Rci is em- 
ployed for the reference length in place of the conventionally used 
nose radius R N . Details of the method of solution are similar to 
those of Refs. 8 and 13 and, therefore, are not presented here. 

Results and Discussion 

Results are presented for the forebody of the SRC shown in Fig. 1 . 
The overshoot entry trajectory 2 (which produces maximum heat 
loads) used in the calculations is given in Fig. 2. The freestream 
conditions at the calculation points for this trajectory are provided 
in Table 2. Peak heating and pressure occur at approximately 54 
and 66 s, respectively, for this trajectory. The SRC forebody is an 
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Table 2 Freestream conditions for Stardust entry trajectory 


Time, 

s 

Altitude, 

km 

Velocity, 

m/s 

Density, 

kg/m 3 

Temperature, 

K 

34.00 

81.64 

12,590.4 

9.63 x 1(T 6 

216.93 

42.00 

71.92 

12,413.4 

1.29 x 10" 5 

221.42 

48.00 

65.44 

12,004.0 

1.06 x lO" 4 

229.00 

54.00 

59.77 

11,136.7 

2.34 x 10~ 4 

238.47 

60.00 

55.02 

9,718.7 

4.39 x lO -4 

248.48 

66.00 

51.19 

7,956.9 

7.21 x 10“ 4 

253.55 

76.00 

46.51 

5,178.9 

1.35 x 10~ 5 

256.90 



Fig. 2 Stardust entry trajectory. 


axisymmetric 60-deg sphere cone with a nose radius R N of ap- 
proximately 23 cm. The computational domain for the present VSL 
results extends to just past the highest point Cl on the shoulder, 
where the flow is predominantly supersonic. The PICA heatshield 
for the forebody extends beyond this point. 2 A 100 x 125 grid is 
employed with all of the computations. Variable grid sizes are used 
both normal to and along the body surface. The minimum distance 
between normal grid points is (2 x 10~ 4 ) R c i . In the direction along 
the surface, the minimum grid size is ( 1 x 1 0“ 1 ) R C \ on the shoulder 
to resolve the flowfield there and is as large as five times this value 
in the nose region to reduce the computational time. These values 
of the grid sizes have been established to ensure grid independence 
of the solution at peak heating condition (r = 54 s). Calculations 
were done by using the Cray C90 computer. The computational 
time required is about 200 and 600 CPU s per global pass for the 
nonequilibrium and equilibrium flow calculations around the body 
(without radiation), respectively. Typically, two global passes are 
required for convergence of the shock shape and surface heating. 

Nonablating Results Along Trajectory with Nonequilibrium 
and Equilibrium Chemistry 

Nonablating stagnation-point total heat transfer rate (conduction 
+ diffusion 4- radiation) is provided at different times along the 
entry trajectory in Fig. 3. Results are obtained with equilibrium 
as well as nonequilibrium (finite rate) flowfield chemistry. Three 
surface recombination boundary conditions, namely, FCW, ECW, 
NCW are used with the finite rate chemistry calculations. Present 
results for the FCW case are in good agreement with those ob- 
tained by Olynick et al. 2 up to time / = 60 s. A maximum value 
of about 1250 W/cm 2 is obtained at / = 54 s from the present cal- 
culations. The corresponding value obtained from Ref. 2 is about 
2% higher. However, the difference between these two values in- 
creases for t > 60 s and may be caused by grid resolution employed 
in Ref. 2. Further, the present finite rate results with an ECW bound- 
ary condition are close to the equilibrium flow results as expected. 
A maximum value of about 1100 W/cm 2 is obtained at t = 54 s in 
this case. The NCW predictions are included for reference and give 
the lowest surface heating with a maximum value of only about 650 
W/cm 2 at t = 54 s. 

The radiation component for total nonequilibrium heating is ob- 
tained from the corresponding equilibrium calculations and is about 
115 W/cm 2 at peak heating (t = 54 s). Figure 4 shows the presently 



Time, s 


Fig. 3 Nonablating stagnation -point total heat transfer rate. 



Fig. 4 Nonablating stagnation -point radiative heat transfer rate. 


computed stagnation-point equilibrium radiative heat transfer rate 
as well as the values obtained in Ref. 2 from a nonequilibrium ra- 
diation calculation. It is not clear why the nonequilibrium radiative 
heating is lower at earlier times and higher at later times in the tra- 
jectory as compared with the equilibrium calculations. Generally, 
the nonequilibrium effects (which are likely to be present at ear- 
lier times in the trajectory) should increase 26 radiative heating in 
comparison with the equilibrium value, and it should approach the 
equilibrium value at later times in the trajectory (with the increasing 
Reynolds number). 

The radiative equilibrium wall temperatures for the heating cal- 
culations of Fig. 3 are shown in Fig. 5. A value of 0.9 is used for char 
emissivity e in the present calculations. Similar to the surface heat 
transfer values, present FCW predictions for surface temperature 
are in good agreement with those of Olynick et at. 2 Reference 2 em- 
ployed a value of 1.0 for e, with zero reflectivity assumed. Figure 5 
also shows Olynick’s values adjusted for e = 0.9. Differences be- 
tween the present predictions and those of Ref. 2 increase at a later 
time in the trajectory for the reasons mentioned earlier. For most of 
the investigated trajectory, the surface temperatures are greater than 
3000 K. Consequently, the FCW boundary condition is physically 
inappropriate because full recombination of air (for FCW bound- 
ary condition) cannot be forced for temperatures greater than about 
2000 K. A physically appropriate surface recombination condition 
for these temperatures is a finite catalytic wall condition, which 
would be bounded by the ECW (most conservative) and the NCW 
boundary conditions. A maximum value of about 3800 K is obtained 
at t = 54 s for the present finite rate results with an ECW condition. 
These results are close to those obtained with the equilibrium flow- 
field chemistry as expected. The NCW calculations give the lowest 
surface temperatures as noted with the surface heating results of 
Fig. 3. 
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£ = 0.9 — o- Equilibrium (Present) 

Finite- rate: 

e = 1.0 — 1 A- Fully catalytic wall (Olynick et al.) 



Fig. 5 Nonablating stagnation -point temperature. 


Ablation Results Along Trajectory with Equilibrium Chemistry 

Figure 6a shows the equilibrium stagnation total heat transfer rate 
with and without ablation along the trajectory. Ablation produces a 
35% reduction in the heating at peak-heating time of t = 54 s. The 
corresponding peak stagnation heat transfer rate of Ref. 2 (without 
ablation) is about 10% higher than the present value, and their re- 
sults also show a reduction of about 35% with ablation. Components 
of the total heat transfer rates of Fig. 6a are shown in Fig. 6b. Re- 
duction in the convective component q c [which is given by Eq. (8) 
and consists of conduction, diffusion, and convection] by ablation 
is caused by injection cooling. Ablation injection reduces the sur- 
face gradients of temperature and that of various species mass frac- 
tions; this causes a decrease in the conductive and diffusive heat 
fluxes. The radiative component, which is relatively small (less than 
11% of the total heating without ablation), is not impacted much 
by ablation injection. There is a slight increase in radiation with 
ablation before the peak heating (/ = 54 s). There is a deeper pene- 
tration of the shock layer by the ablation species C and CO during 
earlier times in the trajectory, and the increase in radiation from 
C line and CO(4 + ) molecular contributions is only partially off- 
set by the absorption of radiation by ablation species during that 
period. 

Surface temperatures used with the equilibrium stagnation heat- 
ing calculations are given in Fig. 7. Surface heating without ablation 
is obtained by using the radiative equilibrium wall temperature as 
already mentioned (see Fig. 5), whereas the ablation temperature 
from Eq. (3) is used for the ablation injection calculations. The no 
ablation temperatures are generally higher (because of the higher 
surface heating) than those with ablation. Present ablation temper- 
ature values are close to those obtained by Olynick et al. 2 up to 
the peak heating time (f < 54 s) in the trajectory. At later times 
the present surface temperature values are lower by a maximum 
of about 18% (at t = 76 s). Also shown in Fig. 7 is the mass frac- 
tion of ablation species at surface, with a maximum value of about 
0.25 at peak heating time (t =54 s). This value implies that 75% 
of the mass at the surface is from the freestream at that time in the 
trajectory. 

The surface ablation injection rate along the trajectory as well as 
the ratio of injection rate to the freestream mass flux, corresponding 
to the heat transfer rate of Fig. 6a, are shown in Fig. 8. A maximum 
value of 3% of the ratio is obtained at time / = 34 s. The value of this 
ratio decreases to about 1.5% at peak heating, where the maximum 
injection rate of about 0.04 kg/m 2 -s is obtained. Even though sim- 
ilar reduction (35%) in heating is obtained with ablation, presently 
computed values of the injection rate and injection mass flux ratio 
at peak heating ( t = 54 s) are about one-half of those obtained in 
Ref. 2. These differences are believed to be caused by the differences 
in the mass fraction of ablation species and their enthalpies used in 
the two calculations. As pointed out in the next section, some of the 
species identified as significant in the present work are not included 
in the equilibrium surface analysis and flowfield chemical kinetics 
model of Ref. 2. 



With ablation 
Without ablatioin 



rate. 

_ Ablation ) 

Radiative equilibrium (no ablation)) Equilibrium *' ow 

— •— Ablation (Olynick et al.); Nonequilibrium How 
C A w : Surface fraction of ablation species; Equilibrium flow 



Fig. 7 Stagnation-point surface temperature and ablator mass frac- 
tion. 
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Present calculations 
Extrapolated 



m 


p v 


Time, s 

Fig. 8 Equilibrium stagnation -point mass injection rate. 



Fig. 9 Stagnation ablator mass fraction and temperature profiles 
[equilibrium flow; peak heating (/ = 54 s)]. 


Laminar/Tbrbulent Ablation Results at Peak Heating 
Condition with Equilibrium Chemistry 

Stagnation Profiles 

Figure 9 shows the temperature and ablator mass-fraction profiles 
at peak heating (/ = 54 s). The effect of ablation injection is limited to 
about 20% of the flowfield close to the surface at that time. This effect 
is also evident from the mass-fraction profiles of the freestream and 
ablation species shown in Fig. 10. The dominant ablation species 
not included in Ref. 2 are C 2 H, C3H, and C 4 H. Noninclusion of 
these species can be partly responsible for the differences in present 
mass loss rate and those of Ref. 2 as mentioned earlier (see Fig. 8). 



Fig. 10a Stagnation freestream and ablation species profiles [equilib- 
rium flow; peak heating ( t - 54 s), /i S h = 0.0105 m]. 


Mole 

Fraction 



Fig. 10b Stagnation ablation species profiles [equilibrium flow; peak 
heating (i t = 54 s), n S h = 0.0105 m]. 


Forebody Surface Distributions 

Surface distributions of the total heat flux for laminar and fully 
turbulent flow cases over the forebody of the Stardust Capsule are 
given in Fig. 1 1. This figure shows the effect of both ablator (PICA) 
mass injection and turbulence on total heating distributions for the 
peak-heating time of t — 54 s along the trajectory. The coupled PICA 
mass injection distributions are shown in Fig. 12. The impact of ab- 
lation injection on total heating is very pronounced for the laminar 
flow. In this case the total heating is reduced by about 35% along 
the forebody (compared with the nonablating surface), essentially 
through the reduction of the convective component as discussed ear- 
lier for the stagnation point. However, for the turbulent solutions, 
where the flow is assumed to undergo instantaneous transition at 
s/R n equal to 0.05, the reduction in total heating is less than 13% 
as compared with the nonablating laminar flow value on both the 
conical flank and shoulder. Obviously the benefit of ablation injec- 
tion in reducing the heating for the laminar flow is partially negated 
when the flow is assumed to be turbulent. The mass injection rate 
distributions of Fig. 12, in general, follow the surface heat-flux dis- 
tributions of Fig. 1 1 . The stagnation nondimensional injection rate 
"i/PooVoo of 0.015 corresponds to a dimensional value of about 
0.040 kg/m 2 -s (as noted earlier). 

There is no noticeable effect of ablation injection and turbulence 
on surface pressure distribution as shown in Fig. 13. The stagnation 
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Fig. 11 Total surface heat-flux distribution [equilibrium flow; peak 
heating ( 1 t = 54 s)]. 
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Fig. 12 Mass injection rate distribution [equilibrium flow; peak heat- 
ing (i t = 54 s)]. 



Fig. 13 Surface pressure distribution (nonablating/ablating) [equilib- 
rium flow (laminar/turbulent); peak heating (/ = 54 s)J. 


(s — 0) pressure in this figure is about 28,000 N/m 2 . This value and 
the pressure distribution are similar to those of Ref. 2. 

Conclusion 

Results are presented for the forebody of the SRC entering the 
Earth’s atmosphere. Solutions are obtained from an axisymmetric 
VSL analysis with and without surface ablation including the effect 
of turbulence. 

The forebody aeroshell consists of a 60-deg sphere cone with a 
shoulder radius one-twelfth that of the nose. For proper resolution of 
the flowfield over the shoulder, the VSL equations are scaled with the 
shoulder radius in place of the conventionally employed nose radius. 
These equations are globally iterated with the Vigneron pressure 
condition to treat the large embedded subsonic region between the 
stagnation line and the supersonic outflow at the top of the shoulder. 

The no-ablation VSL calculations employ an 1 1 -species nonequi- 
librium chemistry model. For these calculations an ECW boundary 
condition is physically consistent and appropriate to use in place of 
the FCW condition (with complete recombination to the freestream 
species). The fully coupled ablation injection calculations are done 
with a 20-species equilibrium chemistry model. With fully equilib- 
rium calculations the elemental conservation equations are solved 
iteratively for each element to determine the elemental composition 
at the surface (adjacent to the ablator) and in the flowfield. In this 
formulation the species boundary condition problem at the surface 
encountered with finite rate calculations is avoided for the case when 
the ablation injection rate is small, and the gas composition at the 
surface is caused by both the freestream and ablation products. The 
small injection rates are usually encountered before large-scale ox- 
idation and sublimation drive the species caused by the freestream 
away from the surface. 

The maximum stagnation heating of about 1250 W/cm 2 is ob- 
tained without ablation injection with nonequilibrium calculations 


and complete surface recombination (i.e., FCW boundary condi- 
tion), whereas a value of about 1100 W/cm 2 is obtained for a more 
realistic ECW boundary condition with a radiative equilibrium wall 
temperature of about 3800 K. Stagnation heating similar to the later 
value is obtained with a fully equilibrium calculation. The maxi- 
mum value of radiative heating component is about 1 1 % at peak 
heating. With ablation injection a decrease of about 35% in the total 
stagnation-point heating (with equilibrium chemistry) is obtained 
at the peak-heating point in the trajectory. Reduction in heating 
is slightly less downstream of the stagnation point and along the 
conical flank, including the shoulder for the laminar case. For the 
turbulent solutions where the flow is assumed to undergo instanta- 
neous transition just downstream of the stagnation line, the heating 
is reduced by only about 13% on the conical flank and shoulder 
as compared with the nonablating laminar flow. Augmentation of 
the convective heating by turbulence appears to partially negate 
the benefit of heating reduction caused by ablation injection in this 
case. 
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